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We construct a two dimensional Cellular Automata based model for the description of pedestrian 
dynamics. Wide range of complicated pattern formation phenomena in pedestrian dynamics are 
described in the model, e.g. lane formation, jams in a counterflow and egress behavior. Mean-field 
solution of the densely populated case and numerical solution of the sparsely populated case are 
provided. This model has the potential to describe more flow phenomena. 



INTRODUCTION 

In pedestrian dynamics, several interesting phenomena 
have been intensively studied: herding effect, lane for- 
mation 1 and jams in a counterflow, egress behavior[2], 
panic situational etc.. Currently two kinds of models 
and their extensions [H[5] have been proposed to describe 
these phenomena: "social force" model pQ and Cellu- 
lar Automata based " floor field" model [7] . The first 
one follows from molecular dynamics model, it involves 
a group of particles interacting by long range "social- 
force" induced by behaviors of the individuals. However, 
the choice of the particular form of long-range "social 
force" is somehow arbitrary and generally pedestrians do 
not decide their behaviour in such a sophisticated way. 
The second one introduced certain kinds of particles that 
transfer inter-person forces. The particle density field is 
often referred to as the "floor field". This model is de- 
rived from models for insects behavior with chemotaxis. 
This model also achieved much success in describing the 
above phenomena, but is somehow a bit far from the re- 
ality. 

In this paper, we present a Cellular Automata based 
model focusing more on how the complex patterns and 
interesting phenomena can emerge from the very simplic- 
ity of the person-person and person-environment interac- 
tion. (From now on, we abbreviate " Cellular Automata" 
to "CA".) This CA model can be regarded as an exten- 
sion of the Ising Model for magnetic spins. Mean-field 
solution for the densely populated case of this model is 
discussed and numerical simulations are carried out for 
the less populated case. Lane formation and jams in a 
counterflow, as well as the egress behaviour are discussed. 

The rest of the paper is organized as follows: In sec- 
tion II the basic CA model setting is presented, a model 
Hamiltonian is given for the system; In section III, the 
extended Glauber dynamics [8] that governs the time evo- 
lution of the system is presented; In section IV, the dense 
populated case is discussed, we apply mean-field approx- 
imation to discuss the equilibrium flow. The results are 
in consistency with numerical results; In section V, the 
sparsely populated case is discussed. We numerically 
simulated lane formation and jams phenomena in coun- 
terflow situation as well as egress phenomenon. Good 



consistency is found with others' results as well as the 
reality. Conclusions are given in section VI. In the ap- 
pendix we present the specific derivation of the mean- 
field solution. 

MODEL 

Our model is based on the local decision making pro- 
cess of a single pedestrian, given its surrounding environ- 
ment and destination direction. The surrounding envi- 
ronment contains the neighbouring pedestrian and their 
walking directions, i.e. pedestrian tend to move in the 
same direction with neighbours; The surrounding envi- 
ronment also contains the available empty spaces: pedes- 
trian tend to avoid crowded area so they move to emptier 
place; they avoid walking towards a wall. The influence 
of the destination direction is modeled by an external 
field that guides them to certain directions. 

We consider a group of pedestrians in certain plain 
space with different destination directions. We call our 
model "N dimensional" if at most a single pedestrian 
can have q = 2N neighbouring pedestrians(see Fig [I]). 
For simplicity we do not consider specific details of a sin- 
gle pedestrian, e.g. his weight, height, age, gender, etc.. 
Since here the space that one pedestrian occupies is the 
smallest space unit we are considering, the space is con- 
sidered to be n discrete lattices, each cell can contain at 
most one pedestrian. Then the total carrying capacity of 
the space is n. The number of pedestrian in this space can 
vary from zero up to n, characterizing the level of crowd- 
edness of the space. In such a N-dimensional model, we 
characterize a pedestrian by q-dimensional vectors 
and ff- q } : s^ 9 * 1 denotes the direction of one pedestrian and 
p( q ) denotes the condition of neighbors of one pedestrian. 
Since our model share similarity with Ising spin model, 
from now on we call a pedestrian a spin. Notice that walls 
and obstacles are equivalent to pedestrians that never 
move. 

The neighbors of one lattice is coded in an order that 
neighbor i and neighbor i + N are on opposite direc- 
tion to this lattice, which then forms a q-dimcnsional 
coordinates. This coordinates system distinguishes with 
N-dimensional conventional coordinates and is designed 
to describe the conditions of the q neighbors of each lat- 
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FIG. 1: schematic drawing of a pedestrian and its surrounding 
lattices of the system, the surrounding lattices are numbered 
so that opposite neighbours are i and i + N,i — N. 



tice. Considering all the possible states of a pedestrian's 
movement, the vectors s^) form a set denoted by S*- 9 -* = 
{eo (9) = (0,0,0,...,0),e! (9) = (1,0,0, . . . ,0), • ■ • ,4° = 
(0,---,l,---,0),---,^ = (0,0,0,. ..,1)}. If g<9) = 
(i — 1, 2, 3 . . . , q), then this spin's direction is towards 
neighbor i, i.e. it will move to neighbor i if neighbor i is 
not occupied by a pedestrian. If = , then this 
lattice is empty. 

Vectors form a set P*- 9 ' = 

{(pi,P2,---,Pi,.--,Pq)\Pi € {1,0}, i = 1,2, 3,. ..,(?}. 
If pj(i = 1,2, 3,..., q) equals one, then neighbor i 
is occupied, otherwise this neighbor is empty. Note 
that only by using q-dimensional coordinates can we 
distinguish between these conditions: 

• Neighbor i is occupied, i.e. Pi = 1, Pi+N = 0. 

• Neighbor i + ^(opposite to neighbor i) is occupied, 
i.e. pi = 0, p i+N = 1. 

• Both neighbor i and i + N are occupied, i.e. pi = 1, 

Pi+N = 1. 

• Both neighbor i and i + N are empty, i.e. pi = 0, 

Pi+N = 0. 

If the total number of pedestrians in the space is M, since 
we apply periodic boundary or the space is large enough, 
then in all the n lattices, the number of occupied lattices 
is M and remains a constant in all process. 

In our model, time is also considered to be discrete. 
It is measured by steps of pedestrians. Again we do not 
consider the details of the difference in footstep length 
and pace of pedestrians. Each step is one time unit, then 
we can introduce the velocity direction vector of one spin 
vi N ) = (si - s 1+N , s 2 - s 2 +n, ■ ■ ■ , s N - S2n) to describe 
the movement of one spin. 

The destination direction is simulated by an exter- 
nal field, denoted by = (B u B 2 , B 3 , . . . , B N ) 
in N-dimensional coordinates. For latter use 
we also introduce a q-dimensional form of 
the external field B^ = (B u B 2 , B 3 , . . . , B q ) 
by defining (B N+i , B N+2 , B N+3 , . . . , B 2 n) = 



— (Bi,B 2 , B 3 , . . . , Bm)- The interaction between 
neighbor spins and interaction between spins and 
external field is described in Hamiltonian of the system: 
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(1) 

The first item results from interaction between the di- 
rections of neighbor spins. Here spins have a tendency 
to be parallel to each other, which represents a kind of 
" cooperation" . If neighbor spins have parallel , then 
they have smaller energy. / > 0, a > are constants 
that determines the properties of the system. This term 
is responsible for herding effect. 

The middle term results from the interaction of one 
spin with its neighbors. Here spins have a tendency to 
move to neighboring lattices which is not occupied. If 
one spin's direction vector points to a neighboring lattice 
which arc not occupied, then v^^ ■ s^^ is 1, otherwise it's 
0. This term results in the tendency to avoid crowded 
situation. The last item results from the external field. 
The spin tends to move in the direction of the external 
field, which simulates its destination direction. 

The Hamiltonian described in the equation can 
be understood directly: herding effect in micro scale 
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is represented by the item —5/^ 
tendency to avoid crowded situation is expressed by 

the item a J2i=i P'^^ an ^ an external field plays 
the role of pedestrian destination direction in the item 
— J27=i B^sf^] The neighboring lattice occupied by 
wall will result in huge increase of energy, then the spin 
will turn away from the wall immediately. 



DYNAMICS RULES 

We do not need to calculate the specific microscopic in- 
teraction force in the dynamical evolution of the system. 
The dynamics here is that the system is approaching a 
energy minimum state. An extended Glauber Dynamics 
[S] is applied, the detailed dynamical evolution rules of 
the CA model is as follows: 

• Denote w(s^ q ^) as the probability per unit time 



of flipping the spin to direction s^ q > 
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E(e^ q ^)is the energy of the system when s^ 9 - 1 — 
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, T is the effective 



"temperature" of the crowd, which reflects the level 
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FIG. 2: schematic drawing of the system, N = 1, 2, 3 



of fluctuation in the system. When temperature 
is high, the spin has higher possibility of chang- 
ing directions, corresponding to an excited crowd 
who are willing to change directions now and then; 
when temperature is low, the spin has lower pos- 
sibility of changing directions, corresponding to an 
inert crowd who tend to keep the same direction. 



ef^ and neigh- 



• If one spin has the direction 
bor i is unoccupied, then move the spin to neighbor 
i. 



• A periodical boundary is applied, i.e. when a spin 
moves to the edge it will appear ont the other side 
of the space. 

We adopt a rule called Q2R Rule [9] developed by Vich- 
niac, basically the rule states that any Cellular Automata 
that changes the state of cells all at the same time can not 
reach equilibrium state, which results from energy con- 
servation law. Under this dynamics, We have the number 
of spins conserved. Numerical realization of the dynam- 
ical evolution share resemblance with Monte-Carlo sim- 
ulation used in various areas of physics as well as math. 
The random number code is from ref|10j. 

In the following sections we will focus on two cases 
of the system. When the lattice is densely populated, 
i.e. vacancy of the lattice only appears occasionally, then 
the homogeneous distribution of spins allow us to apply 
mean-field approximation to probe its dynamical proper- 
ties. When the lattice is sparsely populated, then mean- 
field theory no longer applies, we turn to numerical sim- 
ulation to study counterflow and egress phenomena. We 
show that even though the dynamical rule is simple, the 
phenomena that emerge from the CA model are complex 
and interesting. 



DENSE POPULATED CASE: MEAN-FIELD 
APPROXIMATION SOLUTION 

When the lattice is densely populated, mean-field ap- 
proximation is effective to analyze the equilibrium flow 
of the system(for details see appendix). We consider M 
of the n lattices are occupied by spins, define parameter 
xi to characterize the mean flow, 
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<s N+l - Sl >,l = l,2,...,N 
< si-si- N >,l = N + 1,...,2N 



(3) 



We introduce a mean-field B^f, which contains both 
the influence of the spin-spin interaction and the spin- 
external field interaction. 



B e ffi = fqxi + Bi - a < pi >, I = 1, 2, 3, 



• <1 



(4) 



Then the mean-field free energy of the system can be 
given as: 



F ~ (-fqnj^tf + a < ^ >J2P (q) ~ na< P iq) >) 

(5) 
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Self-consistent Equations 



If the system is in equilibrium, the average of X)"=i sf^ 



can be calculated out: J27=i ^ 
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Also from 
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the definition we have X)"=i ^f' — n < ^ >. If the 
periodic boundary condition is adapted or the system 
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is large enough, then We make a further approximation 
< pi >~< po >! I = 1) 2, 3, . . . , q, consequently 



n < si >= M 



,0(fqxi+Bi) 
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For simplicity and without loss of generality, we only con- 
sider the external field as follows: B\ = B ,Bn + i = 
—Bo,Bk = 0, for other k. Then the solution is 
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Numerical solutions of the self consistent equations are 
in Figure [3^. 

When external field B approaches zero there is still 
non-zero solution when T is below certain critical point; 
Certain low temperature phase transition exists when 
B — » 0+. Approximation of B — > 0+ situation of this 
phase transition is done using Taylor expansion, which is 
showed as the black line in Figure{3a. Define the critical 
point as T c = = l^n' Taylor expansion gives: 
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Since the approximation only holds when temperature is 
close to the critical point, We can see that the black line 
fits well with the numerical solution around the critical 
point and deviates from the numerical solution when T 
continues to decrease. 

Two dimensional simulations are carried out to check 
this mean field approximation results(See Figure |3]b.c.). 
For each parameters we perform 10 separate simulations 
and calculate the average. In order to show the variance 
we plot the maximum and minimum result as well. Gen- 
erally the results fits well, critical temperature turns out 
to be T* ~ 0.4T C . High temperature and low tempera- 
ture limits of the solution is discussed in the appendix. 



Spontaneous Self-organization 
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Simulation results and thoery, 
B=0.025, f=1, a=0, M=7074, n=10000, Tc*~0.43Tc 
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FIG. 3: (a)Numerical solution of xi in Equation Q, N = 1, 
T c = (b)Two dimensional (N = 2) Simulation results 

compared with theory, critical temperature turns out to be 
T c * ~ 0.4T C . (c)Two dimensiona^A^ = 2) Simulation results 
compared with theory, critical temperature turns out to be 
T* ~ 0.4T C . All data are averaged from over 10 separate 
simulations, with maximum and minimum data also plotted. 



Spontaneous self-organization exists when B = 0, then 
the effective field do not contain external field. We re- 
quire the free energy (equation ^) of the system ap- 
proaches the minimum, so =0. If we do Taylor 
expansion on the equation(see appendix), then the ap- 
proximate solution to = is: 

M T I T~ 
x I = ± ] ^^3(l--),i = l I 2,3,... )9 (9) 

The definition of T c is the same with equation ([8]) . The 
uncertainty in the sign as well as in / of xi is caused by 
fluctuations. This is essentially the symmetry breaking 



phenomena in our CA model. The solution is sensitive 
to tiny asymmetry of the initial condition. Physically it 
means: when temperature of the system is below T C) the 
system will form constant flow spontaneously. Note that 
this spontaneous self-organization is different from the 
solution of the self consistent equations when B —> 0+: 
compared with the none external field condition, a tiny 
external field will result in great change in the system. 
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SPARSELY POPULATED CASE: 
COUNTERFLOW AND EGRESS PHENOMENON 

In a sparsely populated case, the total spins in the lat- 
tice M is far smaller than the total number of lattices 
n. Then mean-field approximation does not make sense 
since the system is generally heterogeneous in sparsely 
populated case. So we turn to numerical simulation to 
probe the various phenomena in this case, e.g. lane for- 
mation and jams in counterflow as well as egress phe- 
nomena. 

Counterflow 

Study on counterflow can be found in Refs.pQ (H O \7\. 
Basically, the counterflow situation refers to the en- 
counter of two crowds of pedestrian coming from two op- 
posite directions. The phenomena can take place in sys- 
tem with a periodical boundary or an open boundary and 
also in a narrow corridor. Several distinctive phenomena 
will happen depending on the density of the crowd, the 
average speed of the pedestrian and the crowd's level of 
anxiety. We conducted several simulations of our model 
and showed the different phases of the pedestrian under 
different parameters. Since it's counterflow, in our simu- 
lation pedestrian are labeled into two kinds, either tends 
to be parallel or anti-parallel with the external field. The 
initial positions of pedestrian are randomly chosen. 

Choose the parameters as follows: a — 100, / = 
l,B = l,n = 10000. This corresponds to a 
very active way of walking, pedestrians cooperate and 
also they avoid crowded areas. Figurc[4^ is a snap- 
shot of the initial condition, color indicates the direc- 
tion of pedestrian: black indicates s=(l,0,0,0)(upside), 
green indicates s=(0,0, 1,0) (downside), blue indicates 
s=(0,0,0,l)(right), red indicates s=(l,0,0,0)(left). 

After encounter process(See Figure]!)}) and merge pro- 
cess(See Figure|4p), three kinds of metastable states can 
be observed. We call the three states different phases 
of the system: "block phase" (Figure |4|l), i.e. pedestrian 
of different directions block each other's road and can- 
not move to the other side; "normal phase" (Figure [5^,) , 
i.e. the system is not in order and fluctuates greatly; 
"lane formation" phase(Figure|5^,f), i.e. the pedestrians 
form lanes and move efficiently. Metastable means the 
states keep changing in specific shapes, but remain the 
same kind of pattern. The final state is dominated by the 
total number of spins M and temperature T, the phase 
diagram is shown in Figure]^. Lane formation phase and 
normal phase do not have a clear boundary, lane forma- 
tion phase has a rectangle shape in the diagram approxi- 
mately. A linear boundary between the block phase and 
normal phase is discovered. 

Block phase(Figure |4]i) can also be called jam phase 
in consistence with the terms used by others. In normal 



phase (Figure [5^,) the system is not in order and fluctu- 
ates greatly, though lane formation phase can be seen, 
but the lanes stay just for several time steps and are 
not stable(In each step of time a pedestrian can take no 
more than one step). Transient block phenomenon (See 
Figure [5jD,c,d) is observed before lane formation. In lane 
formation phase(Figure[5^,f), long lanes are formed and 
the system remain relatively stable. Pedestrian can move 
fluently and the transport efficiency is maximum. 

Egress phenomenon 

Egress phenomenon^ refers to the process of a group 
of pedestrian leaving a room in a certain period of time. 
An external field acts as the aim of the crowd that guides 
the crowd to the door. The external field is set so that 
pedestrian are moving towards the up side of the space 
from the exit in the middle(See Figure [6^,). Then the 
crowd clogs around the door until the whole egress is 
over(See Figure [6jD,c,d). During the egress, asymmetric 
behavior of the clogging crowd around the exit is ob- 
served(see Figure [6ji). This phenomenon formed in the 
following process: a pedestrian either from the left side or 
the right side of the clogging crowd enter the exit, then a 
lane is formed since the pedestrian of his side follows his 
movement while the pedestrian from the other side are 
expelled by the lane and await until all the pedestrian 
from the left side enter the exit. 

Figure [7] shows the relationship between egress time 
and door width. Certain critique width is discovered, 
the result is consistent with Ref. j2] , the critical width of 
a door is 3 ~ 4 lattice width. According to Ref. [7], each 
lattice is approximately 40cm x 40cm, which should be 
narrowly able to contain one standing pedestrian, thus 
the critical width of a door is about 1.2m ~ 1.6m, which 
is approximately consistent with our daily life, e.g. typi- 
cal door in Peking University dorm building. 

CONCLUSION 

We construct a Cellular Automata based model to 
describe pedestrian dynamics. We focus on the emer- 
gence of complex phenomena resulting from simple local 
interaction between pedestrian. Analytical mean-field 
solution of the model is given to describe dense popu- 
lated case. Simulations are carried out to verify the re- 
sults. Numerical simulations are carried out to probe the 
sparsely populated case, which is beyond the mean-field 
approximation. Typical phenomena like lane formation 
PQ and jams in a counterflow and egress behavior [2] are 
found. Block phase, normal phase and lane formation 
phase are distinguished and the phase diagram shows 
that only in very limited situation lane formation can be 
observed. A linear boundary between block phase and 



(a) (b) (c) (d) 

FIG. 4: color indicates the direction of pedestrian: black indicates s=(f ,0,0,0)(upside), green indicates s=(0,0,l,0)(downside), 
blue indicates s=(0,0,0,l)(right), red indicates s=(l,0,0,0)(left). (a)Initial state of counterflow, two crowd of pedestrian with 
different directions. (b)Encounter. (c)Merge. (d)Jam or block phase, pedestrian tangle together and cannot form efficient flow. 
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FIG. 5: color indicates the direction of pedestrian: black indicates s=(l,0,0,0)(upside), green indicates s=(0,0,f ,0) (downside), 
blue indicates s=(0,0,0,f )(right), red indicates s=(l,0,0,0)(left). Figure (b)-(f) show the whole process of lane formation. 
(a)Example of normal phase. (b)Lane formation, transient block phenomenon before lane formation happens. (c)Lane forma- 
tion, transient block is starting to approach lane formation phase. (d)Lane formation, transient block is disappearing. (e)Lane 
formation. (f)Lane formation. (g)Phase diagram, To is the numerical unit of temperature 
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FIG. 6: (a)schematic drawing of the external field B and the situation of egress. (b)Initial condition. (c)Egress behavior, (d) 
asymmetric egress, a lane is formed on the left side of the crowd around the door so that the left side go through the door first 
while the right side pedestrian remain blocked. 
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FIG. 7: egress time and the width of the door, critical width 
is between 3 — 4 lattice length 



normal phase is discovered. Asymmetric egress behavior 
is discovered. The critical width for the door is qualita- 
tively measured, results are consistent with Ref.[2] and 
real life. 

Further application of this Cellular Automata based 
model can be expected in qualitative descriptions of other 
flow phenomena. It shows that Cellular Automata is a 
powerful tool to describe non-linear and complex phe- 
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APPENDIX DERIVATION 
Details of Mean-field approximation Solution 

Now we study the hamiltonian of the system 

- n n 



i=l 



we have: 



■JN) -»(j\n / \ i n. 

Vi v j — ( s il — SiN+l,Si2 — Sijy + 2, ■ ■ ■ , SiJV — S i2 N ) 1 ( s jl ~ s jN+l, s j2 ~ s jN+2, ■ ■ ■ , SjN ~ Sj2N) 



Define 



JY 



N 



N 



JY 



— £ Sik ■ Sjk + £ Sik+N ■ Sjk+N — (£ Sik ■ Sjk+N + £ Sik+N ' Sjk) 



k=l 



k=l 



k=l 



k=l 



N 



N 



4 O 4 =; ^2 Sik ■ s Jk+N + ^2 s zk +N ■ sjk 



k=l 



k=l 



It is easy to see that: 



and 

(£q) + #(<?)) Q C {q) = Q C {q) + B {q) Q C (q) 

then 

t N) rt N) = 4 9) -$ ) -4 9) 0$ ) 

Apply mean-field approximation to it, we have 

s f > . s -f = ( S f )_ < s f ) > +s f )) . ( S -f- < s f) > + < g<«) >) 
^ + < 4" ) > - < 4 9) > • < 4 9) > 

since < s| 9 ' >=< >> denote them as < >, then 

sf } • ~ (4 9) + 4 9) ). < ^ > - < ^ > • < S<«) > 

Similarly, 

s f } o 4 9) - (4 9) + 4 9) ) O < s_<<?) > - < > O < ^ > 

^ 9) s f > ~ < S<«> > + < p<9) > . s f )_ < ^) >Q< > 
Then we have the summation as follows: 



*f°^ (JV) - 2 9 < S<«> > £ sf } -nq< S<«) X g<«) > 



<»,j> 



(2g < s< 9) > O 51 4 9) - "9 < > O < >) 



5>1 9) s f > ~< ^ > ]T^ 9) + < > E s 1 9) - n < >< ^ > 

i—1 i—1 i—1 

and it is easy to see 

JV 

< s< 9) > Q < s< 9) > - < s< 9) > • < s< 9) >= - XI ^ lettin S a;/ = (< s; > - < s i+A r >), 



Also we define xn+i = —x\ for later use. 

So the Hamiltonian can be approximated as: 



H -- S effT,$ 9) 

i=l 



1 iV n 

(-fqnJ2 x f + a < >^ / P {q) ~na < p< q) >< s< g) >) 



Here i s tne niean-field, the projection of B (q ^ in direction 1 suffices: 

B effi = fqxi + Bi - a < pi >,l = 1,2,3, . . . ,q 



The partition function of the system is: 



* = E-E 



e -fJH 



~ e 



?<<!) j(9) 
3 1 a n 



since M of the n lattices is occupied by spins, 

Z ~ e -/3(^/9nEr =1 ^+a<^ 9) >Er = iP^ ,) -™a<p<«»>< S -<''»>)^M( e 03 ^^) 

s\ q) eS\{g } 

Free energy 

Free energy of the system is; 

F = ~lnZ 



N n 

{-fqnJ2tf + a< s< 9) > ^ ^ ~ na < ^ >< ^ >) 



J=l i=l 



»feS\{e„} 



Details of Self- consistent Equations 
l.High Temperature Limit Situation 

f3(fqxi + B{] « 1 , so we can do Taylor expansion to the item e^^ qxi+B '^ 

remember x t = (< s t > - < si +N >) = -x N+ iand (B N+1 , B N+2 , B N+3 , . . . ,B 2 n) = -{B 1 ,B 2 

n<st >- M m + Eti ip 2{fqXi + Bir 



so for I = 1 



n < si >~ M 



for other ^ 



2AT + £-=i5/? 2 (/^ + i^) 2 

2^v + E-=i^ 2 (/^ + ^) 2 

1 + /3fqx N+ i 



n< si>~M 



n < s N+ i >~ M 



2A r + ELii/ ?2 (/^ + A) 2 



10 



then 



„ . , M P(fq(x 1 -x N+1 ) + 2B Q ) 

»« * > - < s N+1 » * ^ 2iY + ELi ^ 2(/ ^ + ^ )2 



n« S; > - < Sjv+i » * M m + EUmfqXi + Bi y 
xi = -x N+ i =< si> — < si +N > 

2p(fq Xl + B ) 



nxi ~ M 



2^v + ELi^ 2 (/^ + ^) 2 



2A r + ELi^ 2 (/^ + S,) 2 
then 

fiB 



Xi ~ 



x ; ~0,Z = 2,3,...,iV 



2. Low Temperature Limit Situation 

/3(/«a:, + SO » 1, if (fqxt + B t ) ? 



e 0(fqx,+Bt) 

n < si >= M- 



Guess the solution satisfies 

xi > 0, x; = 0,Z = 2,3, TV 



n < si >= M- 



e f3{fqx l +B l ) 



2N - 2 + eW^i+^o) + e -W9 a!l+ - B °) 
for Z = 1 

e /3(/gxi+B ) 

the solution is 

M 

< si >= — 
n 

Easy to see that for others 

Sl =0, Z = 2,3,...,27V 



so 



M 

.X! = — , xi =0,1 = 2,3,. ..,N 
n 



which satisfies our guess. 
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Details of self- organization 



N 

■3(9) sil) 



2e /3a<P0> . £(1 + + -(/3/^) 4 ) 



2=1 



JV „ N 



~ 2iV e ^<Po> . (1 + £ ±-(0f qXl f + E w(/?/^) 4 ) 



27V 7 ^ 247V v 



ln{ Y e-^'M*)* 



S feS\{e } 



2V N 

ln(2N) + Pa<p > +[£ ^Wfqxtf + E ^(Z^) 4 ] 



2N^ J * ' ^ 247V 
i=i i=i 



N N N N 

~\ E 2lv^^) 2 + E 247V (/3/ ^ )4]2 + 3 E + E ^^^^ 

l=i i=i i=i i=i 



N N N N 

4_2„,2 

2N \hjh-u Z J ^4JV AKn '°" ^^87V 21 """ 
l j=i j=i fe=i 



ln(2N) + I3a<p > +J2^(pf qXl f + ^-l-)(fif qXl f - EE GW*' 1 



N N N N N 

E E 96^/«) 6 ^ + E E E xw 6 ^ 

i=l fe=l 1=1 fclm=l 



F ~ -^n(iC) - j{^(2iV) + Pa < Po >} 



„,.,2 2 N ]\/T N 1 1\/T N N 1 

- ^> E * 2 - f E^x/W + f E E ^/*) V*2 

;=i p ;=i p !=i fe=i 



M N N 1 M N N N 1 

+f EE w«) 6 ^ -f EEE 

P i=l fc=l ' 1 = 1 k=l m=l 

dF_ l Mppg\ M(3 3 fS\ l U 2 , jggVV ^ 2 , MgW ^ 4 

dxf~~ ( 2 IqU 2N > + 47V [ N 3 >Xl+ 47V 2 ^ Xfc + 967V 2 ^ Xk 

l k=l.k=£l k=l 

87V 2 l 6 7V J ^ ^ fc m 

fc— 1 m— 1 



d 2 F _Mp 3 f 4 q 4 1 1 M/3 5 /V 2 
d{x 2 ) 2 ~ 4N ( iV ~ 3> + 48iV 2 Xl 



1 



N 



47V2 l 6 N'^ k 
fc=i 



When /3 > ^ 
From ^ = we have 



so 



or 



d 2 F 



Mf3 z fq A 



dx 2 dx 2 AN 2 



Mf3 5 f e q 6 
48iV 2 



4 



N 



M0*f><f l 1 ^ 2 

6=1 



/I 



_ JV 
3 

1 1 - 



1 



1 
1 

1 1 



N 



1 
1 

1 

1 - 



N 



i 

x 2 



"N-l 



\ *N J 



47V 2 Mf3f 2 g 2 _ fqn 
M(3 3 f 4 q 4 ( 2N 2~> 



( 1 \ 
1 

1 

V 1 J 



Nn 



P 2 pq 2 K Mpfq 



>(1 



-+ / 3 n ^ w ^ 



